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Abstract 



^' Classical sine-Gordon theory on a strip with integrable boundary conditions is 

(— I \ considered analyzing the static (ground state) solutions, their existence, energy and 

stability under small perturbations. The classical analogue of Bethe-Yang quantiza- 
tion conditions for the (linearized) first breather is derived, and the dynamics of the 
ground states is investigated as a function of the volume. The results are shown to be 
C^ \ consistent with the expectations from the quantum theory, as treated in the perturbed 

conformal field theory framework using the truncated conformal space method and 
thermodynamic Bethe Ansatz. The asymptotic form of the finite volume corrections 
to the ground state energies is also derived, which must be regarded as the classical 
limit of some (as yet unknown) Liischer type formula. 

1 Introduction 

Sine-Gordon theory is known as the prime example of integrable field theory, and has been 
the subject of many investigations. In particular, it is well understood on the half line with 
an integrable boundary potential (or boundary condition of the Ghoshal-Zamolodchikov 

type a). 
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In the case of two-dimensional (mostly integrable) theories it is of much interest to 
consider their dynamics in finite spatial volume (strip or circle), since finite size effects (fol- 
lowing the ideas put forward by Liischer [2j for the case of generic quantum field theories) 
provide a very interesting link between different descriptions of the system. Integrability, 
being based on the existence of local conserved quantities, is generally preserved by passing 
from the infinite volume system to a circle, or to a strip with boundary conditions at the 
two ends that are both integrable in a semi-infinite system. Besides that, sine-Gordon 
theory on the strip is also interesting as a description of certain condensed matter systems 

miiniiniizi- 

In the case of sine-Gordon theory with boundaries, such investigations so far mainly 
aimed to connect the exact scattering matrices derived from the bootstrap (infrared data) 
to the perturbed conformal field theory description (ultraviolet data). The main tools used 
to establish this relation are Bethe Ansatz based (TBA, NLIE) or truncated conformal 
space (TCS). While there is a great wealth of information provided by these approaches, 
it is of interest to examine the dynamics of the system from an intuitive, field-theoretic 
point of view, which takes the classical Lagrangian as its starting point and is the goal of 
the present work. 

This approach starts with the classification of classical ground states (vacua) of the 
system, which are identified with static solutions of the field equations with appropriate 
boundary conditions at the ends of the strip. This is followed by an investigation of the 
stability of these states under linearized oscillations of the field. It is very important 
to realize that beyond providing information on stability, these small oscillations can be 
associated with the classical limit of the elementary particle of sine-Gordon theory, which 
is the first breather, therefore their examination gives information on the classical limit of 
one-particle states containing a first breather on the strip. In addition, the stability analysis 
forms the basis of semiclassical quantization of the ground states on the strip (similarly to 
the case of periodic boundary conditions discussed in [H]), but we do not attempt to solve 
this problem in this paper. Instead, we shall be mainly interested in the relation between 
the very intuitive, visual picture of the finite volume dynamics provided by this approach 
to the more abstract, algebraic description provided by Bethe Ansatz and TCS methods. 

The outline of the paper is as follows. In Section 2 we overview some general facts about 
the semi-infinite system, in particular the semiclassical quantization of ground states dis- 
cussed in in, which are needed later. Section 3 outlines general properties of the system 
on a strip, especially concerning the large volume limit. In Section 4 we present a system- 
atic investigation of the case with Dirichlet boundary conditions at both ends of the strip, 
since it allows us to clarify all conceptual issues in a very transparent way. In Section 5 
we present another case with a Dirichlet boundary condition at one end and a special per- 
turbed Neumann boundary condition at the other. Armed with all the lessons we learned 
from these two cases. Section 6 presents a discussion of the generic case. Conclusions are 
presented in Section 7, while two technical issues are relegated to appendices: Appendix A 
contains the description and parameterization of the classical limit of the reflection matrix 
of the first breather, and Appendix B contains some numerical work on the infrared limit 
of the energy levels. 



2 Sine-Gordon theory on the half hne 

The sine-Gordon theory in the bulk is described by the Lagrange density: 

Csg{x, t) = i {d,^x, t)f - V (^(x, t)) 
where the potential is: 

\/ ($) = — (1 - cos/3$) = 2— sin^ ^ > 

Consider the restriction of this theory onto the negative half line, with the (bulk) action 

fO 



/oo /"U 

dt I dx Csci^j't) 
-OD J —OO 



Integrable boundary conditions correspond either to fixing the value of the field at the 
origin (Dirichlet BC) 

$(x,t)U=o = $? 
or introducing a boundary potential of the form (perturbed Neumann BC) 



r{z) 



V^'> {^{z,t)) = M, 



l-cos^($(z,t)-<l>. 



. (2.1) 



The Dirichlet BC with $^ = $o can be recovered from it by taking the Mq — > oo limit. 

The semiclassical analysis of the static solutions in the perturbed Neumann case was 
performed in |2j. Although all the results for the Dirichlet case can be obtained from it in 
the Mq — > oo limit for pedagogical reasons, and because the explicit expressions are needed 
later, we spell out the details here and refer for the general case to ^. 

Using the various bulk symmetries of the model [10] the fundamental range of the 
parameter <l>^ turns out to be 

^ n TT 

< l^o"" = V'o < 2 

Static solutions of the equation of motions satisfying the appropriate boundary conditions 
can be obtained by integrating the static equation of motions: 

Tfl 1 -^ Tfl 

9,'$ = — sin/3$ -^ _(5^$)2 = ^+ (l_cos/3<l>) (2.2) 

Introducing a "time" variable t = —mx and an angular coordinate q = /3$, this system 
describes a pendulum of total energy C. The energy of a static solution is given by 

Emx)] = f dx(^ (9,$)' + V mx))] + Vb ($(0)) 

and the requirement of having finite energy solutions in the field theory fixes the boundary 
condition q{t = oo) = 27m, i.e. C = 0. This system is identical to a classical particle of 
unit mass moving in the potential cosg — 1 as shown on the figure: 




The boundary condition gives the starting position q{t = 0) = 2{pQ. There are two 
solutions satisfying the boundary conditions which describe a particle arriving at one of 
the two nearest potential maxima in infinite time. 

The sine-Gordon analogue of this picture is a soliton (+) or an anti-soliton (-) standing 
at a particular position: 




$+ {x, a+) = ^ arctan L'^i^^-^+A ■ $_ (^x, a') 



27r 4 



arctan e 



where sinh(ma^) = ±cot((/9o). The energies of these two solutions are 

4:771 

E±{ipo) = -^(l=Fcosv9o) 
and they approach their asymptotics at a; — > — oo exponentially fast: 

$+(x,a+) =0(e-™l"l) , <l>_(x,a-) = — + 0(e-™l"l) as x 



-^m{x—a ) 



(2.3) 



-oo. 



(2.4) 



The quantum corrections to these energies can be obtained by analyzing the small 
fluctuations around the appropriate backgrounds: 



(f 



dx^ 



+ m 



2m^ 



cosh (m[x — a^]) _ 



^±{x) = u'^iix) ■ e±(0) = 



The solutions consist of a discrete and a continuous spectrum. 

Introducing e > as o;^ = m^(l — e^) the normalizable modes of the discrete spectrum 
are 

^±{x) oc e™^("-'^*)(e - tanh[m(x - a±)]), 
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where from the boundary condition we have e = =F cos(y9o- Clearly we have a discrete mode 
only for the anti-solitonic ground state which in the light of jH] corresponds to a boundary 
excited state. 

The continuous spectrum can be described in terms of g > (u;^ = m^ + g^) as 

. / N _ ^^^^-iq(x-a±) ^Q + mtanh{m[x - a^]) ^ ^^^iq(^-a±) ^q - mtanh(m[x - a^]) 

iq + m iq — m 

and the ratios A±/B± are determined by the boundary conditions at x = 0. The asymptotic 
form of the fluctuations (x -^ — oo) has the following form: 

e±(x) ^ C±(e'''" + e-*<?"e*^*(«)) (2.5) 

It has an interpretation that the elementary excitation, the linearized form of the first 
breather, reflects back from the boundary with classical reflection coefficient: 

gi5±(q) _ m-iq ±cosypo + f 
m + ig ^ cos (^0 + i 

Observe that the discrete mode can be obtained from the boundary condition dependent 
pole term of the refiection factor and it has an interpretation as the first breather bound 
to the boundary. The other pole singularity is of kinematical origin. 

The semiclassical correction to the energy levels can be obtained by summing up the 
frequencies of these fiuctuations. By performing the appropriate ultraviolet regularization 
[HI at the same time we arrive at the final answer: 

These result, namely the energy of the excited boundary state, the refiection factor of the 
first breather and the energy differences are in complete agreement with the classical limit 
of the quantum results, given in appendix A, if we make the identification ipo = rjd between 
the parameters. 

The system defined on the positive half line can be obtained from the previous one by 
the V : X -^ —X parity transformation. See [TT] for the details of how the two systems are 
related. Clearly the roles of the solitons and anti-solitons are exchanged. 

3 Finite volume: general considerations 
3.1 Static (vacuum) solutions 

Now we consider sine-Gordon theory restricted onto a strip with two integrable boundary 
conditions. The action in the general case has the form 



/oo r nL 

dt / rfx£5G(a;,t) + <^($(0,t))-l4^)($(L,t)) 
oo L^'O 



9$(L,t) 



The last term (proportional to a) is allowed by integrability, and on the strip - in contrast 
to the theory on the half line - it cannot be eliminated from the action by an appropri- 
ate "gauge" transformation [12]. Since this term is a total time derivative, it gives no 
contribution to the classical equations of motion, and we omit it in the sequel. 
The equations of motion and the boundary condition are 



2 

df<^ - dl<!> + -- sin (3<!> = 

(3 

5x$U=. = -^sin|(<l>(z,t)-<l>,) , z = 0,L 

Dirichlet boundary conditions again can be obtained as limits in which M^ — > oo at the 
appropriate end of the strip. 

Ground states (vacua) of the model on a strip are static solutions of the equations of 
motion, i.e. dt^ = 0. It can be easily verified that such solutions provide extrema of the 
static energy functional 

Emx)] = £ dx Q {d.^f + V ($(x))) - <) ($(0)) + vi^) mL)) (3.1) 

Using the first integral ()2.2|1 . the solution, 0(x), can be written in the general form 

X = / , ^ (3.2) 

-^"^'^ ±J2{V{U) + C) 

The sign ambiguity must be fixed from a more detailed analysis which is discussed later. 
The two integration constants $(0) (the value of the field at x = 0) and C are fixed by the 
boundary conditions and the volume. This integral can be written in terms of Weierstrass 
p-function (see explicit examples later). 

It is obvious that £'[$(x)] is bounded from below. In addition, the potential terms are 
also bounded and in fact periodic under $(x) -^ $(a;) +47r//5. These properties guarantee 
that there always exists at least one static solution, but we shall see that in many cases 
there are in fact two different ones (solutions related by a period shift will be considered 
the same). 

Let us consider the limit L ^ oo. In that case finite energy requires that 

$(.) ^ -^ 

for any x far away from the endpoints a; = and x = L. Due to ()2.4|) . we can approximate 
the solution by sewing together two vacuum solutions, corresponding to imposing the 
boundary conditions at x = (x = L) on a semi-infinite system on the half line. In this 
way we get the finite volume solution up to a precision of exp(— mL), (higher precision will 
be described later on), but there is a condition that the asymptotic values of the left and 
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right solutions must match each other. The matching condition leads to some selection 
rules: not all possible pairings of boundary states can be realized on a strip. The selection 
rules can be described as follows. The general boundary excited states can be characterized 
by giving a string of integer numbers 

\ni,n2,...nk) 

(these are the boundary states corresponding to the poles zz/„^, iwn2, ■ ■ -^^ the notation used 
in [nmini)- The two static solutions in infinite volume correspond to the ground state |) 
and first excited state |0). 

Let us first treat the case of Dirichlet boundary conditions. In [12] the authors introduce 
the parity of a boundary state as the parity of the number of integer labels, i.e. of k. For 
the two static solutions, |) has even parity and the asymptotics of the solution is 

$(x) — i> as X -^ — oo (3.3) 

and |0) has odd parity with asymptotics 

<P(a;) —* — as x -^ — oo (3.4) 

This can be extended to the general case since the boundary states \ni,n2, . . .rik) can be 
obtained by binding solitons/antisolitons to the boundary in an alternating sequence. So in 
general even states correspond to configurations with asymptotics ()3.3p . while odd ones to 
(I3.4J1 . Since the configuration on the strip must be continuous, we can infer that for states 
containing no particle the only possible combinations of left and right boundary states are 
either even-even or odd-odd. This is exactly what was seen using truncated conformal 
space approach in |14| . 

If we look at states containing particles, then an obvious extension of the above ar- 
gument yields that states containing even number of solitonic particles obey the same 
matching condition for the boundary states, while those with odd number of solitonic par- 
ticles obey the reverse (i.e. the boundary states at the two ends must have opposite parity 
for these cases). This is meaningful because topological charge is a conserved quantum 
number for Dirichlet boundary conditions. 

The same selection rules can be extended to any boundary conditions, which is simple 
to understand from the fact that the Dirichlet boundary condition can be obtained by 
taking an appropriate limit of the general one. This limit is given in terms of the Ghoshal- 
Zamolodchikov [T] boundary parameters rj and t? by taking ■(9 — *> oo, while the boundary 
state spectrum only depends on rj y^. As a result, it is natural to expect that selection 
rules on the strip must be left unchanged by "undoing" this limiting procedure. The 
rule concerning states containing solitonic excitations can be meaningfully transferred to 
the general case {•& < oo) because topological charge modulo 2 is conserved for general 
boundary conditions as well. 




e +Ri(k)e 



^ik(x-L)^j^^(j^^^-ik(x-L) 



x=0 



x=L 



Ax~l/m 



Ax~l/m 



Figure 3.1: Derivation of the classical Bethe-Yang equation 



3.2 Small oscillations 

We can now investigate small oscillations around a static solution. Putting 

$(x, t) = 0(x) + C,{x) exp{—iu!t) 

(where 0(x) is a static solution) and linearizing in ^(x) we obtain the following equations 
of motion: 






+ V"{cl){x))i{x)=uj^i{x) 



M,^cos^{<P{z)-<^,)i{z) 



, z = 0, L 



(3.5) 



For Dirichlet boundary condition at one (or both) ends we get that ^(x) must vanish at the 
appropriate boundary point (s). Inserting the general static solution (J3.2J1 one gets Lame's 
equation of order 1, with linear boundary conditions for ^(x). 

For the moment let us only consider the limit L —>■ oo. Since in this limit the static 
background is exponentially close to a superposition of two solitons localized at the bound- 
ary, we can apply the asymptotic solution of the semi-infinite case at any point x sufficiently 
far from the boundaries (see Fig. 13.11) . This means that the solution can be written in two 
different asymptotic forms: 

e(x) ~ A (e"'^" + i?i(fc)e''=") as x is large 

e(x) ~ B (e^'^^^-^) + i?2(A;)e-^'=(^-^)) as x - L is large 

where Ri{k) and R2{k) are the appropriate reflection factors calculated on the (semi- 
inflnite) half line ()2.5p . The two asymptotic forms must match in any point x in the 
middle^, which gives two homogeneous linear equations for A and B. The compatibility of 
these equations requires that 



Ri{k)R2{k)e 



2ikL 



^Using standard quantum mechanics, the wave function S_{x) and its derivative S,'{x) must be continuous, 
which translates to the exact matching of the two asymptotic forms for any value of x. 



This is the correct classical limit of the Bethe-Yang equations for a state containing a 
single first breather (i.e. the fundamental particle) in sine-Gordon theory on a strip. This 
general argument shows that the analysis of the small oscillations around a static solutions 
must reproduce the correct semiclassical limit of the state containing a breather between 
the corresponding boundary states. Later we show this explicitly for a number of cases as 
a consistency check of our analysis. 

We remark that one can similarly derive that for periodic boundary condition and a 
background containing a single soliton in finite volume L, the appropriate equation is 

where Ssbi {k) is the classical limit of the phase shift sustained by Bi when scattered by 
a soliton. The derivation in this case proceeds by matching the two asymptotics of the 
scattered wave in a point of the finite volume far away from the location of the soliton. It 
is straightforward to check that this relation is a consequence of the formulae describing 
the small oscillations around a solitonic background in ||Sjj. 

4 Sine-Gordon model on the strip: Dirichlet-Dirichlet 
boundaries 

We start by considering the most simple type of boundary conditions which preserves 
integrability, namely Dirichlet boundary conditions on both end: the value of the field at 
the origin and at L are fixed to be 

$(0) = < , ^L) = ^1 

Using the various bulk symmetries of the model <l>^ can be mapped into the fundamental 
range 

If we would like to describe all the solitonic sectors of this theory we have to consider the 
various boundary conditions 

$(L) = $f + —n ne Z 
simultaneously, thus we will suppose that < $£ < ^. 

4.1 The case $f = 

As a starting point we consider the case when $£ = 0. 

For any n we have one and only one static ground state, which can be obtained by 
integrating the first integral of the static equation of motion (|2.2jl . In the language of the 
classical particle it means that the particle has to reach the top of the hill separating the 
sectors in finite time L as shown on the figure: 

9 




Clearly for this reason in any sector the energy of the classical particle is positive (and 
consequently C > 0). For the sectors 0, —1, . . . the velocity is negative while for 1,2.. .it 



is positive. The actual form of the solution, 0(a;) in the (j){L) = (+) and (j){L) 
sectors are 



27r 



(-) 




mx 



Vo 



/3^/2 



du 



Vsin^ u + C 



mx 



(3^/2 



du 



<P(> 



Vsin^ u + C 



where C 



2m' 



c, 



and C can be determined uniquely from the boundary conditions and the width of the 
strip / as 



/ := m,L 



'PO 



du 



Vsin^ u + C 



(+), I 



du 



v>o 



Vsin^ u + C 



(4.1) 



Note that the solution (+) is in the n = 0, while the (— ) is in the n = 1 topological sector. 
Once C is known, the energies of these solutions can be written as 



4m 



•fa 



E-{vo,l) 



Am 



'sin^ u + Cdu — — 
V sin^ u + Cdu — 



LJipo 



CI 
2 
CI 

T 



(4.2) 



In the large volume limit, / -^ oo, the solution has the form of a standing anti-soliton (+) 
or soliton (-) with the corresponding energies. As < (/9o < f the case (+) forms the 
ground state, while the case (-) is an excited state. These solutions can be obtained by 
parity transformation from the half-line results, since the theory is defined on the positive 
half-line. 



10 



In order to analyze the small fluctuations around this background we rewrite them in 
terms ofWeierstrass's p{z) function. After a straightforward change of integration variables 
and an appropriate deformation of the integration contour into the complex plane the case 
(+) solution can be written as 



dT r dT 

xm 



where 



and 



v/4T3 - g^T - g; J ^AT^-g^T-gs 

Pd Ptb 



Pd = — sm ifo; P^ = — sm — , 

92 = l{l + C + C), g, = -^(1 - C)(l + 2C)(2 + C). 



Since the integrals appearing here are the standard integral representations of p{z) [T 
[TH| we can write 

P^ = ^-^ ~ sm^ ^^ = p{d - mx) , PD = p{d). (4.3) 



dT 
d = .._ ^ (4.4) 



where 

oo 

.= f. 

a/4T3 - g^T - g3 
Pd 

The solution, (/)(x) in the other case (-) is simply 

1-C . 2p<Pix) 

P^ = — sm2 ^1^ = p{d + mx) 

Since C > the roots of the cubic polynomial defining p{z) are 



-C-2 1-C 1 + 2C ,^ , 

es = ^ < 62 = ^— < ei = — ^ — , (4.5) 



thus the half periods of p(2;) are [T^ 



Kik) iK(k') , 1 ,, / C , , 

ViTc Vi + c ViTc Vi + c ^ ^ 

Using the well known relation between p{z) and Jacobi's elliptic functions [121 the classical 
solution </)(a:) can be given explicitly. One finds from ()4.3p in terms of the dimensionless 
xm I— i> X variable 

(4.7) 



/50(x) 
2 


1 f^l + C{dn{dTx,k)±l) 
i V sn{y/l + C{dTx) ,k) 


where d is determined from 






2 1 + C 




"^^ sn^{dy/l + C,k) 
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4.1.1 Infinite volume limit: leading corrections 

It is interesting to consider the / -^ oo limits of the classical ground state energies, 
since they are proportional to the classical limit of the scaling functions ("effective central 
charge"). Using ()4.H) and the asymptotic expansion of elliptic integrals, we obtain the 
asymptotic value of C: 



2 



yC~4 tan^e"' , / -^ oo . 



Plugging this into the formulae for E± ()4.2|] . after some effort we obtain 

E+{^oJ) ~ ^(^l-cos<^o + 4 tan^^e-^') 

^-(¥^0,0 ~ ^(l + cos<^o + 4 cot^^e-^') (4.8) 

It is a very interesting challenge to find a quantum field theoretic derivation of these 
formulae following the general description of finite size effects originally put forward by 
Liischer for periodic boundary conditions [2]. This expectation is supported by the results 
of [H], where, in the periodic case a full agreement is found between Liischer's results and 
the infrared expansion (/ — > oo limit) of the kink's classical energy. 

4.1.2 Stability analysis 

In the semiclassical analysis the field is expanded, $(a;, t) = </)(x) +^(x)e~*'^*, around the 
solution of the static equation of motion, 4>{x). The small oscillations are investigated 
in the approximation when both the equation of motion and the boundary condition are 
linearized. Introducing the dimensionless variable mx i-^ x the boundary conditions and 
the equation of motion ()3.5p take the following form: 

+ 2sin2(M^)e + P^e = O , m = m = 0, (4.9) 

where u^ = m^(l +p^). Using the representation in (I4.3J1 . this can be rewritten as 

^ - [2p{d Tx) + a]ax) = 0, a = -^(1 - C) - p\ (4.10) 

which is the standard form of Lame's equation with n = 1 (cf. |IZ|)- The general solution 
of this equation is known [TH]; if we introduce y = d^ x and ^(x) — ^ z/(y) then it is given 

by 

i/(y) = -^[Aa{y + a)e-^^(") + B(r{y - a)e^^(")], (4.11) 

where A, B are constants, a is the solution of p(a) = a and 

u 

C{u) = / ( p{v) :,]dv, o-(n) =uexp 



U I \ -512 



C(^) --]dv 

V 
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are Weierstrass's ( and a functions. The two boundary conditions on z/ impose the following 
quantization condition on p: 

- 41±lz441±4e--^(") + 1 = (4.12) 

a{d^l + a)a{d- a) 

The algorithm giving the quantized momenta from (j4.12|] is the following: first d and C 
are determined from the boundary parameter $^. (In practice it is simpler to treat C 
as an input parameter and determine d and / from C). Then, using these values in the 



logarithmic version of eqn. (I4.12J1 . one determines a„ for each integer n coming from the 
substitution 1 = exp(2i7m). The n dependent momenta are finally given by 

Pl = -I{l-C)-P{an) (4.13) 

Equation ()4.12|1 (or its logarithmic version) should be considered as the (semi)classical 
version of the Bethe - Yang equation for the first breather (fundamental scalar particle) 
in finite volume l. To support this interpretation we consider the / — > cxd limit, which we 
implement by setting C to zero. Then Ci = 62 = 1/3, 63 = —2/3, and p simplifies [IB] 

2 

p{u) = hcoth^(M). 

o 

Therefore the equations determining a and d become 

cotha = ip , coth^ d = cos^ cpo (4-14) 

Furthermore in this case ( and a can be written explicitly as 

2 
u u 

C(m) = — — + coth(ti), a{u) = sinh(M) exp[— -— ]. 

3 6 



In the asymptotically large / domain one finds from (|4.12|1 

+ 1 = (4.15) 



ip- lip ±cos<^o 2iip 



ip + lip^ cos (fo 

To obtain this we wrote d = do + i^ in accord with (I4.14J1 . where tanh((io) = — cosy^o, since 
the asymptotic background is an antisoliton rather than a soliton standing in the vicinity 
of the X = boundary. In fact we verified numerically that making the substitution 
T = Pd + iy in the integral ()4.4jl yields 

00 



^/4{Pd + lyf - g2{PD + ly) - 93 



do{C,^^)+UJ2iC), 
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where do is real and the purely imaginary uJ2{C) is equal - up to our numerical precision - 
to the second half period of p, which, for C — i> indeed uj2{C) — > if . In the same limit it 
was also found that do —>■ do with taiih{do) = — cosvJo- 

As shown in appendix A the functions multiplying e^*'^ in eqn. ()4.15p are the (semi) classical 
limits of the first breather's reflection amplitudes on various Dirichlet boundaries with pa- 
rameter (fio in infinite half space. This confirms the interpretation of eqn. ()4.12j) as the 
classical limit of the finite width Bethe-Yang equation. 

The essential difference between the two cases in eqn. ()4.12[ I4.15|] is that for the lower 
sign - that corresponds to the background describing the excited ground state - the function 
multiplying e^'''^^"'^ (e^*'^) admits a pole, that for / — > cxd may be interpreted as a bound 
state, while no such pole is present for the upper sign - that corresponds to the ground 
state background. 

In the half line case poles of the reflection functions provide information on the discrete 
spectrum, which corresponds to excited boundary states. Now the potential poles come 
from the quantization conditions ()4.12p . 

Certain poles of the functions multiplying e^'^*^"-* in the various quantization conditions 
may correspond to bound states of the flnite width problem. Other poles of the same 
functions may describe some kinematical singularities. The difference between these two 
types of poles is in the way they depend on the boundary parameters and on the width of 
the strip: the potential bound state poles depend on both the boundary parameters and 
on /, while the kinematical poles are independent of them. Furthermore the (semi)classical 
energy of the bound state in flnite strip 



/ 2(1 — C) 
uj = mA/l+p2 = mdl ^— — --pia) (4.16) 

should go, for large /-s, to the classical limit of the energy of the known bound state in the 
half line theory. 

Consider the quantization condition (I4.12J1 . With appropriate signs in the arguments 
the function appearing here has a potential bound state pole and the energy of the corre- 
sponding bound state is 



tu^ = myC + sin>o- (4.17) 

Note that for C ^ (/ ^ oo) u^ -^ msimpo indeed. The bound state is present only for 
the excited, solitonic, background, which is consistent with the / -^ oo asymptotic. 

4.2 The general Dirichlet case 

Now we relax the condition ^fj* = and describe the general case. Keeping all the possible 
sectors of the theory we allow the ^{L) = $f + ^ boundary conditions. For any n we 
have a sector of the theory. The static ground states can be visualized in the language of 
the classical particle as shown on the flgure 



14 




7^=^ 




For the sectors labeled by —1, —2, . . . there exists a unique static solution. The energy 
of the particle is positive while the velocity is negative. For the sector 1,2,... we have the 
same statement: the static solution is unique, the energy of the particle is also positive 
just as is its velocity. The n = sector is the most complicated one. Here the energy is 
negative but we do not have a unique static solution. The two lowest lying static solutions 
correspond to the following motions of the classical particle: 




In the large volume limit / -^ oo the energy of the particle tends to zero from below 
and the solutions have interpretations as two independent soliton and anti-soliton standing 
at the two ends of the strip, i.e. we recover the solutions corresponding to two separate 
Dirichlet boundary conditions corresponding to $^ and ^^. In the case (A) the solution 
has asymptotic energy E+^ipo) + E^{ipi) while for case (B) it is E_{(^q) + E_{{pi) (see 
eqn. ()2.3|) for notations). Since < y^o < | the case (A) solution has lower energy if 
(po < 71 — ipi. We can restrict ourselves to this case without loss of generality, since the 
other possibility can be described by parity transforming the strip and correspondingly 
inverting the time variable of the particle, which exchanges the roles of the initial and 
final conditions. (Notice, however, that increasing continously ipi over it — ipo the roles 
of the two solutions, namely the ground state and the first excited states are exchanged). 
Now decreasing the volume the particle has less and less energy; when the kinetic energy 
becomes exactly zero at the one of the endpoints of the motion (ipo or (pi) then the (B) 
type solution ceases to exist and no such solution exists for smaller volumes. At that point 
a breather comes out from the boundaries and this time-dependent solution will follow 
the energy line as it will be explained later. We can conclude that although there are 
more static solutions in this sector for general volume, only the one with the lowest energy 
survives in the small volume limit. 
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The mechanism described above is clearly analogous to what was seen before in [T^fTO] . 
where in a different configuration an excited boundary state could "decay" in finite volume 
by emitting a breather. The appropriate scaling function was described by the Bethe-Yang 
equation of a single-particle state containing a first breather with zero momentum quan- 
tum number below the critical volume, and by continuing the rapidity of this particle to 
imaginary values above the critical volume. The phenomenon described above is somewhat 
different, since in this case the two boundaries can be thought of as emitting a soliton and 
an antisoliton, respectively, which then combine to form the breather inside the strip and 
it is not obvious how to describe it in terms of some Bethe-Yang equation as was done in 



4.3 n = topological sector: non monotonic background 

4.3.1 (A) type solution 

Now we consider the SG model on a strip with Dirichlet boundary conditions at both ends 
with: 



0<(^o < 



TT 



, < v?; < VT - vJo, 



(4.18) 



and consider the case (A) solution first. In this case, for sufEciently large /-s, the ground 
state is given qualitatively by the superposition of a static anti-soliton and soliton located 
in the vicinity of the two boundaries: 



2n 



x=0 



x=L 



Therefore the solution first decreases from $^ to some minimal value e then increases 
from e to $P and in these two domains it is given by slightly different expressions. Since 
dx^ must vanish at the turning point the integration constant C is 



2,/5e. 



C7 = -sin^(^)<0. 



Using the dimensionless x variable the solution can be written as 



vo 



X 



dv 



(3<P/2 



Vi 



sm^ V 



C 



e < (/) < $ 



D 
' 



16 



in the decreasing domain, while as 

X = l- + 



dv 



Pe/2 



Vsin^ V + C 



e < 6 < ^ 



D 



in the increasing one. Here we introduced 



dv 



fie/2 



Vsin^ V + C 



and the turning point is determined by the width of the strip as 

l = l_ + L = 



dv 



dv 



Vsin^ V + C J Vsin^ v + C 

/3e/2 f3e/2 

Note that I — > oo corresponds to e — i> 0. The classical energy of the solution can be written 
as 



E 



Am 



V 



- sm^ — + 
2 2 



/ du^/sm'^u + C + / du\^sm'^u + C 



lie/2 



fie/2 



I 



The large / limit of this energy has a simple interpretation as the sum of the energies of 
the static solitons/antisolitons constituting the ground state. 

In both domains the solution can be expressed in terms of the Weierstrass's p(-2). The 
difference now is that since C < the roots of the cubic polynomial defining ^{z) are 



es 



-C - 2 1 + 2^ \-C 

< 62 = — - — < ei = 



3 ~ "'^ 3 ~ "' 3 ' 

thus the half periods oi^{z) are [12] 

uJx = K{k), uj2=iK{k'), k = Vl + C, k' = V^. 

In deforming the various integrals we also encounter the following expression 



(4.19) 



(4.20) 



dT 



Pe 



l-C 



sm 



./3e. 



1 + 2C 



Since ^^^ = 62, and p{uji + UJ2) = 62 thus I = cui + uj2. Therefore, defining di i = 0,1 as 

previously 

1 — C 
Pd, = — ^ sin^ ipi = p{di), 
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one can write 

I = do + di- 2{uJi + 1^2), I- = do- {uji + UJ2), 1+ = di - {ui + UJ2), 



and the solution becomes 
in the decreasing domain and 



P<j> = Pido -x), 
P^ = p{x -do + 2{uJi + UJ2)) 



in the increasing one. 

Using the relation between p{z) and Jacobi's elliptic functions [11] the classical solu- 
tion 0(x) can be given explicitly in this case also. One finds from (14. 3|) in terms of the 
dimensionless xm 1-^ x variable 

/30(x) _ 1 j^ n±cn{d-x,k) \ 



2 i \ sii((i — X, k) 

where d is determined from 

cos^ (po 



sn'^{d, k) 

In the stability analysis the fluctuation equations are solved simultaneously in the two 
domains by exploiting the p{u) = p{—u) and p(u + 2{uJi + 002)) = p{u) invariances. Using 
the coordinate y = do — x the boundary condition at x = becomes 



u{do) = Afido, a) + Bfido, -a) = 0, /(x, a) = ^i^±^e-^«") 

a[x) 



•) 



while at X = / one finds 

iy{-di + 2{uji + U2)) = Af{-di + 2{uji + UJ2), a) + Bf{-di + 2{uJi + cua), -a) = 0. 

Non vanishing A, B exist only if the 

2lC(a) ^i(^i ~ " ~ 2(^1 + Ld2))(r{do - a) ^ ^ 
a{di + a - 2{iUi + LU2))(r{do + a) 

condition holds. Using the quasi-periodicity of a{u) and an interesting identity [11] one 
can write 

a{di - a - 2{uJi + uj2))a{do - a) a{di - a) 

■exp [AaC[uJi + 002)] 



a{di + a — 2{iUi + uj2))crido + a) a{di + a) 

a'^{a + ui + uj2)cr{di — a) 



(4.22) 



a'^{-a + uji + uj2)(T{di + a) 

and this converts the quantization condition into 

^2ic(a) ^ (^0 -/- + «) tT {di -l+ + a)a [dp - a) a{di - a) ^ ^ 
a (do — l_ ~ a) a {di — l^ — a) a (d^ + a) a{di + a) 
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The advantage of this form is that the / -^ oo limit (when also l± -^ oo) can be carried 
out easily using the asymptotic forms of the Weierstrass functions described earlier. 

As stressed earlier the solution described so far is valid for sufficiently large I. As I 
decreases from infinity e increases from zero, and the qualitative nature of the solution 
changes when e reaches the smaller one of the two boundary values ^^in- Indeed /+ 
becomes zero when e = $^jjj and this particular /q is the smallest strip width when the 
ground state solution has an extremum; for / < Iq the solution becomes monotonic and 
there is only one domain. In this single domain the solution can be written as 



Vo 



X 



dv 



P4>/2 



/3«!>/2 



Vsin^ V + C 



, or X 



dv 



Vsin^ V + C 



90 



where C > — sin^ ipQ and is determined by the width of the strip 



v'o 



dv 



fi 



Vsin^ V + C 



or / 



dv 



Vsin^ V + C 



fi 



fo 



When expressed in terms of the Weierstrass's p function this solution again takes the 
Pij) = p{do — x) form but / is expressed now as I = do — di. The stability analysis in the 
single domain leads in a straightforward manner to the quantization condition 



cr{do + a)a{di — a) 



I </o. 



4,3.2 (B) type solution 

Next we consider the SG model on a strip with the same Dirichlet boundary conditions as 
before 

, <(/?;< TT 



< v?o < 2 



V'o, 



but look for the (B) type solution. In this case, for sufficiently large /-s, the solution first 
increases from <l>^ to some maximal value x = ^ — e then decreases from x to ^f and in 
these two domains it is given by slightly different expressions. 




x=0 



x=L 
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Since dx^ must vanish at the turning point the integration constant C is 

C = -sin^(^) = -sin^(f)<0. 

Repeating the same steps as in the case of the (A) type solution one easily obtains that 
this solution is expressed as 

P^ = pido + x), P^ = p{x + do + 2{uJi + U2)) 

in the increasing (resp. decreasing) domain. The quantization condition for the fluctuations 
around this background is obtained by the di — > —di i = 0, I substitution in eqn. (|4.23l) . 



This, in the classical limit, leads to a change of sign for the cos ipi terms. 

The solution described here exists for sufficiently large L As / decreases e increases and 
the existence of the solution comes to an end when e reaches the smaller one of the two 
boundary values $min- ^^ contrast to the (A) type solution at this particular / = /min both 
/+ and /_ remain finite and the solution does not exist for / < /mm- 

To understand this consider the poles of the functions multiplying e^'^*^"-* in the quan- 
tization condition as functions of /. For this excited ground state there are two bound 
state poles, which depend on both the boundary parameters and on /; the energies of the 
corresponding bound states are given by ()4.17|) with the (po ^ ipi i = 0, I substitution. It 
is clear that one of the cj^'s becomes zero at the minimal volume, i.e. when C becomes 
C = — sin^ v^min- The appearance of this zero mode clearly indicates the instability of the 
underlying excited ground state solution and is in accord with the picture described at the 
end of Section 4.2. 

4.4 The —1 topological sector: monotonic background 

Next we consider the sine-Gordon model on a strip with DD boundary conditions charac- 
terized by the same (po and (fi boundary parameters as before. The fundamental range of 
the boundary parameters is as in ()4.18p . which describes the zero charge sector. In the —1 
topological sector, the right boundary value of the field can obtained from the value in the 
zero charge sector by the following shift: 

271 



D 



$(L) = $^ 



P 



In this case, for sufficiently large / the static ground state is the superposition of two 
anti-solitons rather than the superposition of a soliton and antisoliton. 
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Because of the different signs, in the effective point particle description the particle has 
to go over the maximum of — V^($) at $ = 0. As a consequence the integration constant 
C has to be positive C > and the solution becomes monotonic. 

In this case there is no need to introduce l± and the two domains but to emphasize the 
similarities between this case and the previous one we define 

dv r , f dv , , X , , , 

, .. : = uJi + uj2-do; l+= / , .- z = di-{uji + uj2); I = l_+l+. 

Vsin^ v + C J V sin^ v + C 



Using arguments similar to the case in subsection 4.1.1, we obtain the following result for 
the large volume asymptotics of the energy (i.e. the classical limit of the scaling function): 

E (v9o, (fij) ~ -^ (2- cos (po + cos (fi + S tan — ^ tan — - — ^e"' j (4.24) 

It is easy to see that the leading correction vanishes ii ipo = or ipi = it , which is consistent 
with (pTSll . 

In both domains the ground state solution is given hj P^j, = p{x + do). Using this in 
the stability analysis leads to the quantization condition 

cr((io — a)a{di + a) 

It is difficult to compute the / -^ oo limit, as the background does not show explicitly the 
presence of both anti-solitons. To circumvent this define do instead of do: 

/_ = / , .^ = do- {^1 + ^2) =uji + uj2- do; 

J Vsin^v + C 




then, using the identity ()4.22p . the quantization condition can be recast as 

2ic(a) ^(do -1- + a)(j{di -1+ + a)a{do - a)(y{di - a) _ 
(7{do — /_ — a)(r{di — /+ — a)a{do + a)a{di + a) 

The C ^ 0+ limit, when a\\ I, l± ^ 00 can now be computed easily. 

4.5 UV limit of ground state scaling function 

The / — > limit of ground state energy can be calculated simply by observing that in small 
volume the energy functional ()3.ip is dominated by the kinetic term 

^ 1 



2 



d.M 
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so the ground state solution (with Dirichlet boundary conditions) is simply 
The energy of this solution is given by 

^ (*f - *?)' 

which agrees with the result from conformal field theory (see formula C.7. in [TO]). 

The fact that the ground state solution is unique in the UV limit is consistent with our 
earlier observation that one of two static solutions that exist for large volume is destabilized 
as the volume decreases. 

5 Sine-Gordon model on the strip: mixed boundaries 

Now we consider a mixed three parameter family of integrable boundary conditions on the 
strip: we assume Dirichlet boundary condition 

$(0) = $0^ 

at the lower end of the strip, while at x = L the perturbed Neumann boundary condition. 
Using the bulk symmetries of the theory $^ can be mapped into the fundamental range, 
while $L can be in the domain: < $l < % and Ml > 0. 

In this section we consider a simpler form of the general boundary condition (SGB): 

a,$|i = <5Mi^sin(^), 5 = ±1, (5.1) 



depending only on the M^ parameter, which is assumed to be non negative, and the positive 
(negative) sign in front corresponds to choosing <l>^ = ^ ($^ = 0) respectively. Below 
we sometimes consider an even simpler case when ipo = /?$^/2 = 0, called Neumann type 
(NT) problem. 

5.1 Classical static solutions 

At X = L we may express C in terms of ^l = ^{L) using the boundary condition: 

C = -{l-A ]——sm^— = ——C, A =— . 

^ ^ P^ 2 P^ Am 

Note that C is negative for A^^ < 1, i.e. as M^ grows slowly from zero but stays in the 
[0, Am/ 0^) domain. As {dx^Y must be non negative, this may impose some restrictions on 
the classical solution. In particular at the x = Dirichlet end 

sinVo-(l-^"')sin2^>0. 
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Thus, for A^^ < 1 in the NT problem {ipo = 0) $l must stay either zero or 27r//? indepen- 
dently of L, which implies that the solution must be the constant $(x) = 0. (Note that 
for A~^ > 1 this conclusion does not hold). 

The classical solution (j){x) is obtained from ()3.2I1 by carefully correlating the sign in 
the boundary condition with that originating in the square root: 

mx= -^ 5 = -l(A), mx= -^ , 5 = 1 (B) . (5.2) 

i/3#/2 Rw y^o R[u) 



Here R{u) is the positive root 



/3$i 



R{u) = W sin^ M - (1 - A-^) sin^ '—^ = Vsin^ u + C, 

and (f){x) < ^Q in ()5.2I A), while 0(x) > <l>^ in ()5.2I B)^. The equation determining 
(f){L) = ^L in terms of the dimensionless length of the strip / = mL has the form 

W /3*i,/2 

f du , , ^ f du , ^ 

/3*l/2 (^0 

in the two cases. In case (A) / changes from zero to infinity as ^l decreases from $^ to 0. 
The same happens to / in case (B) as $l increases from $^ to 27r//5, if sin^ ipQ> 1 — A^^. 
If, however, sin^ ipo < 1 — A~^, then, in case (B), ^l must be bigger than $niin (/5'^min > tt), 
which is obtained from 

sm v9o = (1 — ^ )sm — - — , (5.3) 

and consequently the solution exists only for /min < / < oo. 

The energy of these classical solutions can be written in all cases as 



/3'J'l/2 



£=|^ (l_^-=)!,ta^^+ , I fl(„)rf„ + ,^-.eos^ 



(5.4) 



fo 



Interestingly, the second term is positive, while the last term, which is the contribution of 
the boundary potential, is negative for large l-s for both the (A) and the (B) cases. For 
asymptotically large l-s the energy has a simple interpretation as the sum of the energies 
of the static (anti)solitons building up the ground state. Indeed for the case (A) say it is 

A-7T) 

E^ = -^ {--^'^ + 1 - cosv^o) , 

where the first term is the energy of a special perturbed Neumann soliton in the left half 
space while the second and third terms give the energy of a Dirichlet antisoliton in the 

^In these considerations we always assume that < $^ < 7r//3. 
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right half space. The asymptotic expression of the energy in the (B) case differs from this 
only in the sign of the cosine term in accord with the fact that now there is a Dirichlet 
soliton in the vicinity of the boundary at x = 0. 

The energy of these classical solutions as functions of the dimensionless width of the 
strip I was investigated numerically for a wide choice of parameters ^q and A. These 
investigations showed that these E{l)-s are smooth functions that tend rather rapidly to 
their asymptotic values given above; for 6 = —1 (case A) this happened already at / ~ 2 
while for (5 = 1 (case B) at / ~ 8 — 10. Furthermore these investigations revealed that 
if for a given / both solutions exist then always (i.e. not only asymptotically) the energy 
of the (A) type solution is the smaller. This fact is in accord with the idea that the (A) 
type solution is the classical ground state and the (B) type one is some sort of classical 
excitation. 

Particularly interesting is to compare the energy E{1) of the (B) type solution of the 
NT problem (<l>^ = 0), when A~^ > 1 and that of the constant $ = solution. Actually 
even in this case the (B) type solution exists only for / > Zmin where 



/•^*^/2 du , lA-^ + 1 , , 

lynin = lim / -^^-^ = In \ —— (5.5) 

At I = /min the first two terms vanish in the energy expression (15. 4|) thus E (Zmin) = ^A~^ 
which is precisely the energy of the constant solution. The numerical study showed that for 
^ > ^min E{1) is always smaller than this particular value (see Fig. 15. H where the continuous 
line denotes E{1), x denotes its asymptotic value and the o stands for the energy of the 
constant solution). 

We can rewrite the solutions in terms of Weierstrass's p{z) function as (I4.3J1 with d 
given by ()4.4|) . For < A^^ < 1 the parameter C < so the roots of the cubic polynomial 
defining p{z) are given by 1)4.19^ and the half periods by ()4.2flj) . In the case when A~^ > 1 
(C > 0) the corresponding expressions are given by ()4.5|1 and ()4.6p . Using the relation 
between p{z) and Jacobi's elliptic functions |Tn| the classical solution 0(x) can be given 



explicitly. When C < and ^-^ = ei one has ()4.21L while for C > and ^^ = 62 we have 
(I4.7J1 . Finally we mention that in case of the (B) type solutions one has to change d — x to 
d + X in eqns. ()4.3I |4.21l , I4.7|l to get the representation in terms of elliptic functions. 

5.2 Stability analysis of the NT problem 

Before considering the more general cases we briefiy describe the stability analysis of the 
0(a;) = static solution. This solution exists for all values of ^~^, when A~^ < 1 this is 
the only solution, when A~^ > 1 and 6 = 1 then there is also an ()5.2I B) type solution. 

In the stability analysis we write $(x, t) = + ^(x)e~*'^* and linearize both the equation 
of motion 
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Figure 5.1: E{1) in the NT problem with A'^ = a/S/^ 



and the boundary conditions 

m 



0, 



eix)i 



s^^m 



4 



around = 0. In case of the discrete spectrum we write cu^ 
solution of the linearized equation of motion is 



m 



then the general 



From the two boundary conditions it is found that A, B can be different from zero only if 

e-(5^"4anhe/ = 0, I = mL . (5.6) 

For 5 = —1 this equation has only the e = solution and this implies vanishing ^, thus 
no instability of the classical solution. For 5 = 1 this equation has a nontrivial solution if 
I > ^crit = ^; for / — > oo this solution tends to A~^: e -^ A^^. If A^^ < 1 this implies 
the existence of a bound state, while for A^^ > 1, when u^ becomes negative at large /, 
it means an instability. In the latter case the instability starts at /q where e(/o) = 1, from 
this condition it is easy to show that Iq = l^am (cf. eqn. ()5.5|l ). This observation clarifies 
the meaning of the instability at I > Iq: it signals the transition from the constant solution 
to the (B) type one which has a lower energy. 

In the semiclassical quantization (201 the energy of the boundary bound state associ- 
ated to the nontrivial solution of eqn. ()5.6|) is given hy lu = m^X — €"^{1). Particularly 
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interesting is the / dependence of this bound state energy, for large / one finds from ()5.6|1 
(0 < A-^ < 1) 



2A 



~2 



cu(Z) = m(Vl - ^-2 ^ — _^ exp(-2^~^/)), for / ^ oo. 

This shows that the finite width of the strip gives an exponential correction to the bound 
state energy in the infinite system. Based on our experience with semiclassical quantization 
of boundary sine-Gordon theory in infinite half space we expect [H] that the (semi)classical 
limit of an exact, fully quantum computation of the finite size correction of the boundary 
bound state energy would coincide with the expression above. 

In case of the continuous spectrum we write oj"^ = m^(l + g^) which converts the solution 
of the equation of motion into the sum of plane waves. The two boundary conditions now 
impose the following 'quantization condition' on q: 

oi„iq + i6A^^ 

1 + e -^ ^nr^ = 0- 

q - t6A-^ 

There are several things that should be remarked about this equation. First of all the 
expression multiplying the exponential is nothing but the product of the classical limits 
of the first breather's (B^) refiection factors on a ground state Dirichlet boundary with 
$^ = (which is trivial = 1) and on an SGB boundary. In fact S = —1 corresponds to B^ 
refiecting on the ground state boundary |), while 6 = 1 and A~^ < 1 to B^ refiecting on 
the first excited boundary |0) (see appendix A). 

The second remark about the quantization condition concerns its relation to the finite 
size correction to energy of the boundary bound state. For this we recall that in infinite 
half space the bound state shows up as a pole in the reflection amplitude at a purely 
imaginary value of the rapidity and use the quantization condition to determine the shift 
of this pole for large but finite l-s. To this end we put q = sinh6' and continue 9 to 
pure imaginary values 9 = iv. Furthermore, for A^^ < 1, introducing sinu = A~^ the 
quantization condition with 6=1 becomes 

sin w + sin u 

1 + exp(— 2/ smf J = U. 

sin V — sin u 

For large / from this equation we find 

2 sinu 
V = u — exp(— 2/sin'u) h C'(exp(— 4/)), 

COSM 

and this implies that the leading correction to the bound state energy is 

2sin^M / ^, . N 
mcosv = mcosu + m exp(— 2/ sm-u). 

COSM 

This of course coincides with the previous result, but now it is possible to interpret the 
coefficient of / in the exponent and the coefficient of the exponential correction in terms 
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of the classical limit of the reflection factor. In fact in these simple considerations the 
classical nature was never exploited, thus had we started with a (quantum) Bethe-Yang 
equation for a single particle with a (quantum) reflection amplitude having a bound state 
pole at V = u with residue iG we would have obtained for the bound state energy at finite 
/ 

mcosv = mcosu — msinu G exp{—2l sinu) . (5-7) 

We think that this expression is completely general and would give the leading finite size 
correction to the energy of any boundary bound state in the half line theory. We also 
expect that a derivation of this formula must exist based on Liischer's ideas [2]. 

5.3 Stability analysis 

In the general case we write $(a;, t) = 0(a;) + ^(x)e~"^*, where (f){x) is the solution in ()5.2I 
A,B). Introducing the dimensionless variable mx i— > x the equation of motion take the form 
as in (14. 9|] but with the boundary conditions: 

e(0) = 0, e'(a:)b = 5^-^cos(^^^e(0, 

For the (A) type solutions, using the representation in ()4.3p . this can be rewritten as ()4.in|l 
for which the general solution is given by ()4.1HI . The two boundary conditions on u impose 
the following quantization condition on p: 

a{d - I - a)a{d + a) _^i^^^-^ -C{d ~ I) + ({d - I - a) + ({a) + Q ._^ .. „^^ 
aid-l + a)a{d - af -({d - I) + ({d - I + a) - ({a) + Q ^ ' ' 

where 

n A/1-1 f^^L 
Q = oA cos ' 



=-1 

The quantization condition for the fluctuations around the (B) type solutions is obtained 
by making the d -^ —d and 6 = 1 substitutions in eqn. ()5.8|] . 

To support the classical Bethe-Yang interpretation of eqn. ()5.8p we consider the / — ^ oo 
limit, which is implemented by setting sin^ ^y^ (and consequently also C) to zero and 
perform the same calculation we did in the Dirichlet case. Since ior I ^ oo Q —^ — -^^, in 
the asymptotically large I domain one finds from (15. 8j) 

_ ^p - Up + cos ^,^,,^^p + A-^ ^ ^ ^ ^ 
ip + lip — cos (^0 —ip + A ^ 

The asymptotic quantization condition for the fiuctuations around the (B) type solutions 
is obtained by changing the signs of cos^Jq and A~^ in ()5.9jl : 

ip + lip — cos (fo —ip + A~^ 
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The important difference between eqn. ()5.9|1 and ()5.in|) is tfiat tlie expression multiplying 
the exponential in latter one admits such poles that, for / — * oo may be interpreted as 
bound states while no such pole is present in (15.911 . The classical energy of these bound 
states is given by 

uji = msiny^o, ^2 = mVl — A'"^. 

It is shown in appendix A that the expressions multiplying e*^'^ in eqns. (|5.9|5.inp are 
the (semi)classical limits of the first breather's reflection amplitudes on Dirichlet and SGB 
boundaries in infinite half space. 

For the (B) type solution there are potential bound state poles in the quantization 
condition when < A~^ < 1 coming from both the Dirichlet and the SGB ends. The 
energy of the bound state associated to the Dirichlet boundary is given by eqn. ()4.17p . 
while the energy of the bound state belonging to the SGB boundary ()5.H1 is given by 



uo 



N ,..,„„„ /^-^ 



m\ cos \\/\ ~ A" 



In this case again oj^ = is obtained at the critical value of / (/min), when the static 
solution in question ceases to exist since then ^l becomes ^ram thus C + siv? lpq = — (1 — 
A^"^) siv? ^ ^'° + sin^ (yJo = in view of eqn. ()5.3|) . The interpretation of this zero mode is 
the same as in the DD case (cf. end of Section 4.2). 

6 General boundary conditions 

Now we discuss the ground state and fluctuations in boundary sine-Gordon model when 
on one of the boundaries the most general, two parametric boundary condition, while on 
the other the simplest Dirichlet b.c. is imposed: 

9.$(a;)|o = -Mo^sin {^ ($(0) - O^ ; ^{L) = 0. 

The integration constant appearing in 

can be expressed in terms of the parameters characterizing the boundary condition at x = 



as 



C = (a~' sin (f ($(0) -K))- sin ^) (^"^ sin (^ ($(0) - K)) + 



sm 

2 y 



Note that C > is necessary to guarantee that $'|l is real; we assume this in what follows. 
Repeating the steps of the previous derivations it is easy to see that the dimensionless 



28 



width of the strip determines now the x parameter as ^ 

/3<J.(0)/2 

du 

t 







Vsin^ u + C 



The advantage of this form of C is that it makes obvious the asymptotic ( / -^ cxd) limit 
of $(0): for this one must have C —> thus 



2 A-^cosS^±l 

In the following it is assumed that $^ > and $(0) > thus $(0)^^ is used as the 
asymptotic limit. 

The classical ground state 0(x) can be written now as 

1-C . 2/50(x) , , ^ , , 1-C . 2/5*(0) 
— sm — - — = p[n — X); where p[n) = — sm — - — , 

and the stability analysis leads in a straightforward way to the quantization condition 

a{n - I + a)a{n - a) _2;c(«) Q " C(^) + C(^ - «) + C(«) ^ ^ 
a{n - I - a)a{n + a) Q - (^(^n) + ({n + a) - ({a) 

where Q = ^^^ cos | ("^(0) — $^). The large / limit of this quantization condition is 
obtained by replacing the Weierstrass functions with their C ^ limits; this way, after a 
long but straightforward algebra that involves some surprising cancellations, one finds 



^2ilp!P_ 



1 A-^ cos ^ - p2 + ipJi + ^-2 + 2^-1 cos ^^ 



^P + ^ A-^ cos ^ - J92 - ip^l + A-^ + 2^-1 cos ^ 

Comparing with P| one can check that this is consistent with the classical reflection factor 
on a general integrable boundary. 

7 Conclusions 

The aim of this paper is to provide some intuitive picture of boundary sine-Gordon theory 
on a strip by investigating the classical dynamics of the theory. First of all, it is observed 
that a matching rule exists in flnite volume, meaning that not all combinations of boundary 
states are allowed on a flnite strip. We have seen that this can be extended to all boundary 
states (not just the ground states described by static solutions) and the resulting rule is 
perfectly consistent with that observed in TCS approach. 

^This form is valid if x > 0; for x < one should interchange the upper and lower limits of the integral. 
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It was also established that in general for large values of the volume there are two 
different classical ground states, but under certain conditions one of them is destabilized 
as the volume is decreased. We discussed the nature of this instability and it turned out 
that it corresponds to the appearance of a breather (of zero momentum quantum number), 
formed by a soliton/antisoliton emitted by the left/right boundary (or vice versa). We 
have shown that this behaviour perfectly agrees with the TCS spectra, and also with the 
quantum Bethe-Yang equations. In particular, for three cases we derived the asymptotic 
form of the finite size correction for the energy of the ground state ()4.8|4.24p . and we also 
gave the leading correction for a boundary excited state in (I5.7J1 . It is an interesting open 
question to derive these (and similar formulae) from an extension of Liischer's argument 
(originally formulated for periodic boundary conditions in [2]) to the strip. Work is in 
progress in this direction. 

The results of the present work lay the foundations for a full semiclassical quantization 
of ground states on the strip, which we leave for further investigations. Such results could 
be used to test and complete NLIE descriptions proposed for the description of finite size 
effects, both in the case of Dirichlet [IHIEII and more general boundary conditions fT2| . 
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Note added. 

After the submission of this work a paper appeared [25| containing results that partially 
overlap with ours. 

A The classical limit of B^ reflection amplitudes 

In the general case the reflection factor of the flrst breather, -B^, on the ground state 
boundary |) is given by [2T] 



B\i\9) = ''! ^'' ',, ^ )\' . ^i)l^ . V. , (^) = .. )l -Li ■ (A.l) 



,(1)^.. _ (I) (^ + ^){^.~ \) (^ - \) ... _ -nh (f + .f ) 

Here 6 is the rapidity of 5^, A = ^ — 1, while rj, -d are the two parameters characterizing 



Stt 
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the bootstrap solution p. (5^'s reflection factor on the flrst excited state |0), R.J (9)^ is 
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obtained from this expression by the substitution r] ^ r] = 7r(A + 1) — r] [IHl)- In [HI it is 
shown that in the (semi) classical limit A — i> oo one should scale the r] and -d parameters 

77 = r/ei(l + A), ^ = ^,1(1 + A), 

and keep rjd and i?ci finite. Furthermore it is also shown, that as long as < t^cI < f the 
only potential pole in the ground state reflection amplitude is out of the physical strip thus 
B^ cannot create a bound state on |). On the other hand B^ can create a bound state on 
|0) since in the same semiclassical limit RU{9) has a pole at ~ f — ?7ci, which is in the 
physical strip. 

This analysis applies to both the Dirichlet case - when the reflection factor is obtained 
in the ■(9 — *> 00 limit of ()A.Hl - and the simpler general (SGB) one ()5.1j) - when the 
reflection factor is obtained either by the ■(? = {A~^ < 1) or by the r] = {A~^ > 
1) substitution from (|A.Hl . Using the relation between the bootstrap and Lagrangian 
parameters [H US! 122], one finds in the Dirichlet case 

thus as long as < $^ < | there is no bound state in the ground state reflection while 
there is one in the reflection on |0). From the general expressions derived in P| it follows 
that in the SGB case 

cosr^ci = A'^ if A~^ < 1, and cosh'i^ci = A~^ if A'^ > 1. 

Finally starting from ()A.1|1 and using p = sinh 9 instead of the rapidity one finds in the 
(semi) classical limit 

DT}i^)f^-, _ ^p-1 ip + COS r]d 
ip + 1 —ip + cos ?7ci 

in the Dirichlet case while 



'B\^\p) = z , : V, :^ (A.2) 



^^(1) JP + cosr^ or ^4\p) = ^P±^^^ (A.3) 

'' ~ip + cosr/ci '' —ip + cosh-yci 

in case of the SGB boundary. The expressions for the limiting values of the reflection 
amplitudes on the flrst excited state are obtained by making the substitution rjd -^ tt — 77^1, 
which amounts to changing the signs of the cosine terms in eqns. (jA.2|A.3|l . This completes 
the interpretation of the expressions appearing in the asymptotic quantization conditions 

B Numerical analysis of the quantization condition 

In this appendix we investigate numerically the difference between the predictions of the 
exact flnite volume quantization condition and the asymptotic (/ -^ 00) one to get an 

31 



estimate on the accuracy of the latter. In the analysis we consider the DD boundary 
conditions with ipi = and the (+) case only. 

In the procedure (with a fixed v^o) the input parameter is C, from this / and d are 
determined first, according to ()4.HI and ()4.4|] . then, using them in the logarithmic version 
of (IKl^ 

2/C(an) + log — 7- ^7;^ — -, r]=2z7in, B.l 

an is obtained, and finally the square of quantized momentum is found from (I4.13JI . The 
asymptotic quantization condition is obtained from ()B.lj) by replacing the Weierstrass 
functions with their C ^ (I ^ cxd) limits described earlier; this way one obtains 



2il%)n + log 



{%Vn- 1) (ipn + COSygp) 

(ip„ + 1) {i'Pn -cosv^o) 



2mn . 



(B.2) 



We carried out the numerical analysis with y^o 
in the table below: 



7r/6; the results for n = 1 are summarized 



c 


I 


v' 


pIs 


{j^llf 


0.4 


0.757175 


17.079765 


16.876796 


YlflX^'^ll 


0.1 


1.287219 


5.843204 


5.774261 


5.956549 


0.01 


2.376994 


1.680518 


1.671835 


1.746801 


0.001 


3.523447 


0.760057 


0.759118 


0.794994 


0.00081 


3.628742 


0.716545 


0.715781 


0.749527 



Here p^ is the momentum squared obtained from eqn. (IB.2|1 and the last column shows 
the momentum squared of a free particle. These result indicate that the predictions of the 
asymptotic quantization condition are very good estimates of the exact finite volume ones 
even at intermediate strip widths. (Note that here / is measured in units of the fundamental 
scalar particle/first breather mass / = mL while in TCSA studies it is usually measured in 
units of the soliton mass ^XCSA ~ ^^ ~ [Ml"^^"^^ — {M/m)l). 
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